The origin of the positron excess in cosmic rays 
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We show that the positron excess measured by the PAMELA experiment in the region between 
10 and 100 GeV may well be a natural consequence of the standard scenario for the origin of 
Galactic cosmic rays. The 'excess' arises because of positrons created as secondary products of 
hadronic interactions inside the sources, but the crucial physical ingredient which leads to a natural 
explanation of the positron flux is the fact that the secondary production takes place in the same 
region where cosmic rays are being accelerated. Therefore secondary positrons (and electrons) 
participate in the acceleration process and turn out to have a very flat spectrum, which is responsible, 
after propagation in the Galaxy, for the observed positron 'excess'. This effect cannot be avoided 
though its strength depends on the values of the environmental parameters during the late stages 
of evolution of supernova remnants. 



The PAMELA satellite began its three-year mission 
in June of 2006, and among its goals there was that 
of measuring the spectra of cosmic ray positrons up to 
270 GeV and electrons up to 2 TeV,each with unprece- 
dented precision [l[ . Recent results [2] show that the ra- 
tio of positrons to electrons plus positrons (the so-called 
positron fraction) in the cosmic ray spectrum appears 
to rise with energy, at least up to ~ 100 GeV, as al- 
ready found by previous experiments, including HEAT 
[1| and AMS-01 0], although with smaller statistical sig- 
nificance. The clear discrepancy between the observed 
positron fraction and the predictions of the standard 
model for the origin and propagation of cosmic rays in 
the Galaxy, led to many possible explanations, ranging 
from the annihilation of non-baryonic dark matter [5| to 
the possibility that new astrophysical sources, especially 
pulsars [6j , could provide the additional positron flux. It 
is worth recalling that both lines of thought lead to the 
"correct" spectral slope rather naturally, but they are 
very different in terms of providing the correct normal- 
ization of the positron fluxes. The dark matter interpre- 
tation typically requires large, and somewhat artificial, 
annihilation rates. Such large rates could, in principle, 
result from dark matter possessing an annihilation cross 
section in excess of the value predicted for a simple s-wave 
thermal relic [av «3x 10 -26 cm 3 / sec), for example due 
to the Sommerfeld effect Q , or from non thermal WIMPs 
[H . In the case of pulsars on the other hand, the energetic 
requirements appear to be all but extreme, although an 
efficiency factor needs to be introduced by hand [6j, and 
the mechanisms for escape of the pairs from the pulsar 
environment are basically unknown. 

At this point it is useful to recall what are the pre- 
dictions for the positron flux of the so-called 'standard 
model' of the cosmic ray propagation in the Galaxy: cos- 
mic rays are assumed to be accelerated in astrophysi- 
cal sources, such as supernova remnants (SNRs), with a 
source spectrum Ncr{E) cx E~ 7 (throughout the paper 
we adopt the convention that for any particle distribution 
function / one has 4irp 2 f(p)dp = f(E)dE). The power 
law behaviour, in the case of SNRs, naturally arises from 
diffusive shock acceleration 10]. The spectrum of cos- 



mic rays observed at Earth is ucr{E) cx Ncn(E)T esc (E), 
where r esc (E) is the escape time from the Galaxy. This 
time scale is inferred from the ratio of Boron to Carbon 
(B/C) fluxes in cosmic rays and is typically found to scale 
as T eS c(E) cx E~ s , with S ~ 0.3 — 0.6. In this standard 
model, positrons only arise as secondary products of cos- 
mic ray interactions in the Galaxy. Positrons are there- 
fore injected at a rate Q+(E) cx ncR(E)nnca (where 
a is the cross section for the relevant process), while 
the equilibrium spectrum is n+(E) cx Q + (E)r e (E) cx 
E^- 5 T e {E), where T e {E) « Min[r esc (E), n oss (E)], and 
ti oss is the loss time scale during propagation in the 
Galaxy (this scaling may however be affected by the spa- 
tial distribution of the sources with respect to the diffu- 
sion region). In the range of energies we are interested 
in, in most cases energy losses are dominant over escape, 
therefore r e (E) rs ti oss (E) cx 1/E (losses are dominated 
by inverse Compton scattering (ICS) and synchrotron 
emission). It follows that the spectrum of positrons ob- 
served at Earth is n+ (E) cx E^- 5 - 1 . The injection 
spectrum of electrons in SNRs is usually parametrized as 
N e (E) = K ep NcR(E), so that following the same rea- 
soning the equilibrium spectrum of primary cosmic ray 
electrons is n_(E) cx E^^ 1 . It is straightforward to see 
that the ratio n + /n^ cx E~ s , namely it decreases with 
energy, at odds with PAMELA findings. These simple 
arguments are nicely summarized in [9(. In this paper 
we show that the standard model as summarized above 
lacks one important phenomenon which is intrinsic in the 
acceleration process. We also show how taking into ac- 
count this phenomenon a positron excess arises in the 
standard picture of acceleration of cosmic rays in Galac- 
tic sources. The spectrum and spatial distribution of cos- 
mic rays in the accelerator are obtained by solving the 
transport equation [Tol . [Tl| and one finds the well known 
result that the solution is 



fcR,(x,p) = K ( — 
\Po 



F(x,p), 



(1) 



where p is the particle momentum, F(x,p) = 1 down- 
stream (x > 0) and F(x,p) = exp(ux/D(p)) upstream 
(x < 0). The predicted cosmic ray spectrum at the shock 
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surface (x — 0) is a power law in momentum and the 
slope is 7 = iui/ (ui—u 2 ) = 3r/(r— 1), where u 12 are the 
fluid velocities upstream and downstream and r — u\ju 2 
is the compression factor. For a strong shock r — > 4 and 
7 — ► 4. For relativistic particles (p ~ E) the spectrum 
in energy is Ncr(E) = <iirp 2 fcn(p)u 2 T S N ~ £ ,_7+2 , and 
for a strong shock Nqr{E) ~ _E~ 2 , the well known result. 
Here rs at is a typical age of a SNR. 

Cosmic rays accelerated at the shock produce sec- 
ondary + e~ inside the source through hadronic in- 
teractions with production and decay of charged pions, 
while the decay of neutral pions leads to production of 
gamma rays. It is crucial to realize that this process 
occurs in the same spatial region around the shock in 
which cosmic rays are being accelerated. Therefore it is 
unavoidable that secondary e~ take part in the same 
acceleration process. Notice that these particles are al- 
ready suprathermal, therefore we do not need to worry 
about their injection. The production rate at a position 
x around the shock is 

Q±(x,E)= J dE'N CR (E' ,x) d<7{ f E E) n aas {x)c, (2) 

where c is the speed of light, n gas is the gas density for 
pp scattering in the shock region (n gaSi2 — rn gaSi i) and 

da± dE' ^ 1S ^ ne differential cross section for a proton of 
energy E' to produce an e + or an e~ of energy E. Here 
we calculate these cross sections following the approach 

of [3. 

The calculations are carried out neglecting the non- 
linear effects which lead to the appearance of a precursor 
upstream of the shock (see [l3[ for a review) . Moreover 
we do not comment here on some important aspects of 
the escape of particles from the remnant which are as 
relevant here as for ordinary cosmic rays, but imply a 
difficult integration over time of the acceleration history. 

The most important point of this paper is the calcula- 
tion of the equilibrium spectrum of e~ + e~ produced in 
the acceleration region. Their behavior is described by 
the transport equation which automatically takes into ac- 
count the presence of the shock, the advection with the 
fluid and diffusion: 

df± n , ,d 2 f± ldu df± 

U ~d^ = D{p) ^ + zTx p -dp- + 0±(x ' p) ' (3) 

where f±(x,p) is the equilibrium distribution function 
of electrons (-) and positrons (+) (number of particles 
per unit volume per unit energy). In the equation above 
we neglect the effect of energy losses of electrons because 
here we will be restricting ourselves to situations in which 
the secondary pairs have energy below that for which 
synchrotron and ICS losses are important. We will check 
a posteriori that this is the case. 

Eq. ([3]) has to be solved with the following bound- 
ary conditions: f±(x,p) and df±/dx vanish at upstream 
infinity (a; — > — oo) and df±/dx remains finite at down- 
stream infinity (x — > +oo). The boundary condition at 



the shock surface is obtained in a straightforward way by 
simply integrating Eq. ^ in a narrow region around the 
shock. With these conditions, one can easily show that 
the solution to Eq. © is in the form 

f±(x,p) = f±,o(p) + — x, (4) 
u 2 

where /±,o(p) = f±( x = 0,p) and is the solution of the 
following equation: 

P-Q^ = -lf±,o + l(-^+r )—^-Qi(p), (5) 

where D\ (p) oc p a is the diffusion coefficient upstream of 
the shock, and Qi(p) — Q±(x = 0~ ,p) is the rate of in- 
jection of pairs immediately upstream of the shock. The 
factor £ ~ 0.05 represents the mean fraction of the en- 
ergy of an accelerated proton carried away by a positron 
or electron in each scattering. 

The solution of this equation is promptly found to be: 




Qx{p'). (6) 



The physical meaning of Eq. ((4]) is that secondary parti- 
cles which are produced within a distance ~ D{p)/u from 
the shock (on both sides) participate in the acceleration 
process (first term in Eq. (HI)), while the downstream 
secondary particles produced farther away are simply ad- 
vected and their density increases with x (second term 
in Eq. ^) up to a maximum distance ~ u 2 tsn, where 
tsn is the age of the SNR. The contribution of SNRs to 
the flux of secondary e~ and e + is the integral of Eq. (|4]) 
over < x < u 2 tsn- 

One should keep in mind that roughly Q\{p) oc 
fcn(p) ~ (at least at p > lGeV/c), therefore Eq. 
dU leads to /±, ~ P~ 7+Q , where a > is the slope 
of the diffusion coefficient as a function of energy: the 
equilibrium spectrum of the pairs that take part in the 
acceleration is flatter than the injection spectrum of sec- 
ondary pairs. This is crucial to understand that the flux 
of secondaries in the sources cannot be simply rescaled 
with grammage from the interstellar medium cosmic ray 
interactions because of the strong dependence of positron 
production on the spectrum of primaries. 

At this point we can calculate the equilibrium spec- 
tra of primary electrons and secondary e" + e~ in the 
Galaxy after propagation. We model the Galaxy as a 
cylinder with half-height H = 3 kpc. We adopt a diffu- 
sion coefficient in the Galaxy D ga i(E) = D (E / GeV) s , 
with S = 0.6 and Dq chosen so to have an escape time 
at 10 GeV of T esc (10GeV) « 20 Myr: D Q rj 10 28 cm 2 s- x . 
Energy losses of electrons and positrons are due to ICS 
and synchrotron emission and we adopt a typical total 
rate of energy losses b(E) = w 2 x \Q- 1& E 2 GeV GeV/s. 
The simplified picture that arises is leaky-box-like with 
all limitations that it implies. 

If TZs/S is the rate of supernovae occurring in our 
Galaxy per unit volume, the equilibrium distribution of 
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cosmic rays in the Galaxy is 



ncB,(E) = N C RllsN"r eac (E). 



(7) 



The equilibrium spectrum of secondary e~ + e~ pro- 
duced by cosmic ray interactions in the Galaxy is de- 
termined by a balance between injection, losses and 
escape from the Galaxy. For the diffusion coefficient 
D{E) ss 10 28 Eq® v cm 2 s _1 the loss time is shorter than 
the escape time at all energies above ~ 10 GeV, namely 
at all energies of interest for us. In this case the equilib- 
rium spectrum of the diffuse secondary pairs can easily 
be written as 



n±(E) 



K N n H c 
b(E) 



dE h 



,w da ± (E',E") 
dE n C R{E )- 



(8) 



dE' 



where n# is the gas density averaged over the volume 
of the Galaxy (including disc and halo) and a coefficient 
.Kjv ~ 1-2 — 1-8 is introduced to account for the inter- 
action of nuclei other than hydrogen. Following [T3| we 
use Kn = 1.8. Clearly, the choice of a different diffu- 
sion coefficient in the Galaxy may lead to the need for a 
more detailed solution, taking into account the interplay 
between escape and losses. Moreover if a non-leaky box 
model is used, a slightly different slope of the equilibrium 
spectra is obtained, though the positron fraction remains 
unaffected. 

Similarly, for the secondary pairs produced inside the 
sources, one has: 



i a ± {E) = K n Hsm 



b(E) 



dE'N s ± (E'), (9) 



where N s ± (E)dE = 4tt P 2 [/±, + (1/2)Q 2 t S n} u 2 T SN dp is 
the distribution function of pairs at the sources in energy 
space instead of momentum space (we integrated Eq. §Z§ 
over the downstream volume, exactly as for CRs). 

Finally, for the spectrum of primary electrons in the 
sources we adopt the standard procedure of assuming 
that N e (E) = K ep N C R.[E), where K ep « 7 x 1CT 3 . The 
equilibrium spectrum of primary electrons is then: 



,{E) = K ep K S M 



b(E) 



dE'N CR (E'). (10) 



Before illustrating the results of our calculations we dis- 
cuss briefly the choice of diffusion coefficient in the accel- 
erator, which is not the same as in the Galaxy, because of 
the generation (and damping) of turbulence in the shock 
region, either due to the same accelerated particles [llj 
or due to fluid instabilities. Here we carry out the cal- 
culations for a Bohm-likc diffusion coefficient, which we 
write as: 

D B (E) = K B -r L (E)c = Z.Zxltf 2 K B B- l E Ge v cm 2 s -\ 

(11) 

Here B M is the local ordered magnetic field in units of 
\iG and the coefficient K B — (B/SB) 2 allows to consider 
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FIG. 1: Positron fraction as a function of energy. The data 
points are the results of the PAMELA measurement. 



faster diffusion (K B > 1), which is common when mag- 
netic field amplification is not as efficient. 

These are all the ingredients needed for the calcula- 
tion of the positron and electron fluxes at Earth. The 
positron fraction, defined as the ratio of the total flux 
of positrons to the total flux of e~ + e + , is plotted in 
Fig. [TJ The data points are the results of the PAMELA 
measurement. The error bar on energy is of the order 
of half the distance between two consecutive data points. 
The solid line refers to the case of maximum energy of 
the accelerated particles (and therefore also of the sec- 
ondary particles after reacceleration) E max = 100 TeV, 
while the dash-dotted and dotted lines refer respectively 
to E max = 10 TeV and E max = 3 TeV. The dashed curve 
represents the standard contribution to the positron frac- 
tion from secondary diffuse pairs. We adopt a reference 
age tsn ~ 10 4 years for a SNR. The three curves refer 
to {K B ,n gas ,uB^u s } = {20, 1.3, 1,0.5} for E max = 100 
TeV, {20, 2, 1, 0.5} for E max = 10 TeV, and {20, 3, 1, 0.5} 
for E max — 3 TeV (n gaSy i is the gas density close to the 
SNR in units of 1cm -3 and u% = ui/10 8 cm/ s). One can 
see that these values are appropriate for old supernova 
remnants, which however are also expected to be the ones 
that contribute the most to the cosmic ray flux below 
the knee. Unfortunately during such phase the maxi- 
mum energy of accelerated particles decreases in time in 
a way which is very uncertain: slowly in the case of no 
damping and rather fast if effective magnetic field am- 
plification and damping are present. This is the reason 
why in Fig. [1] we considered the three values of E max . 
A solid evaluation of this effect can only be achieved by 
carrying out a fully time dependent calculation (Caprioli 
and Blasi, in preparation). A prediction of this scenario 
is that the positron fraction grows and eventually levels 
out at ~ 40 — 50%. The fluxes of electrons and positrons 
are plotted in Fig. [2] for the case E max = 100 TeV. We 
assumed that the closest source of cosmic rays is located 
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FIG. 2: Fluxes of e" and e + at Earth for E max = 100 TeV. 
The dotted line refers to primary electrons, the dashed lines 
are the fluxes of positrons (upper curve) and electrons (lower 
curve) from interactions of cosmic rays in the Galaxy. The 
dot-dashed lines are the fluxes of positrons (upper curve) and 
electrons (lower curve) from production in the sources. The 
thick solid line is the total flux. The data points are from 
Fermi/LAT QJ]. 



at a distance of order ~ 1 — 2 kpc, so to introduce a high 
energy cutoff at ~ 1 TeV, namely when the propagation 
time from the closest source exceeds the loss time (this is 
a strong function of the distance to the closest source). 
A cutoff may also be produced by the acceleration pro- 
cess in the sources. Recent observations by ATIC [15j . 
HESS [IB] and Fermi/LAT seem to confirm the presence 
of the high energy cutoff in the diffuse electron spectrum. 
The dotted line refers to primary electrons, the dashed 
lines are the fluxes of positrons (upper curve) and elec- 
trons (lower curve) from interactions of cosmic rays in the 
Galaxy. The dot-dashed lines are the fluxes of positrons 
(upper curve) and electrons (lower curve) from produc- 
tion in the sources. The data points are from Fermi/LAT 
(Tc| . A few remarks are in order: 1) at energies above 



~ 20 GeV the main contribution to the positron flux is 
strikingly the one of secondary pairs in the sources. 2) 
The flat spectrum of the secondary pairs, energized inside 
the acceleration region, makes them provide up to 50% 
of the total flux of electrons (and positrons) at Earth at 
high energy. Their contribution is sufficient to flatten the 
total e + + e~ spectrum, which is the quantity effectively 
measured by ATIC and Fermi/LAT. 3) No sharp feature 
appears in the total spectrum, at odds with ATIC data, 
which arc however not confirmed by Fermi/LAT results. 

We also stress that a contribution to this positron flux 
might come from a fraction of SNRs located in proxim- 
ity of dense molecular clouds, where the target for pp 
collisions may be enhanced and the importance of the 
mechanism made more evident. These SNRs might also 
have sufficiently high surface brightness (despite the old 
age) to be detected by Fermi/LAT, provided the target 
density for pp interactions is large enough. 

We can conclude that the positron excess can be a 
consequence of acceleration of cosmic rays in SNRs or 
other sources. The scenario discussed here has numer- 
ous implications: first, a flux of antiprotons is predicted 
to be produced, which is compatible with present data 
pjj; second, a similar effect, though more model depen- 
dent, appears in the spectra of secondary nuclei which 
are now being calculated. As far as individual sources 
are concerned, one has to keep in mind that no effect 
is expected from young, bright SNRs. In old remnants a 
possible signature might appear in radio spectra but may 
be easily confused with other effects. In any case super- 
novae of type la are more likely to contribute to this effect 
because they typically explode in denser media. 



Acknowledgments 

The author is grateful to R. Aloisio, E. Amato, D. 
Caprioli, G. Morlino and P. Serpico for endless discus- 
sions and exciting collaboration. 



[1] |htt p: / / pamela.roma2. Infn.it / index. php| 
[2] O. Adriani et al. Preprint larXiv:0810.4995l 
[3] S. W. Barwick et al. [HEAT Collaboration], Astrophys. 
J. 482, L191 (1997); S. Coutu et al. [HEAT-pbar Collab- 
oration], in Proceedings of 27th ICRC (2001). 
[4] M. Aguilar et al. [AMS-01 Collaboration], Phys. Lett. 

B646, 145 (2007). 

[5] I. Cholis et al, Preprint larXiv:0809.1683l 
[6] D. Hooper, P. Blasi and PD. Serpico, JCAP 1, 25 (2009); 
I. Biisching, I., O.C. de Jager, M.S. Potgieter and C. 

Venter, Astroph. J. Lett. 678, 39 (2 008). 

[7] M. Cirelli and A. Strumia, Preprint larXiv:0808.3867l 
[8] P. Grajek, G. Kane, D. J. Phalen, A. Pierce and S. Wat- 
son, Preprint larXiv:0807. 15081 



[9] PD. Serpico, Phys. Rev. D79, 1302 (2009). 
[10] R. Blandford and D. Eichler, Phys. Rep. 154, 1 (1987). 
[11] A.R. Bell, MNRAS 182, 147 (1978). 
[12] T. Kamae et al, Astroph. J. 647, 692 (2006) [Erratum: 

Astroph. J. 662, 779 (2007)]. 
[13] M.A. Malkov and L. O'C. Drury, Rep. Prog. Phys. 64, 

429 (2001). 

[14] M. Mori, Preprint larXiv:0903.3 260 

[15] J. Chang et al., Nature (London) 456, 362 (2008). 

[16] F. Aharonian et al. [HESS Coll.], Phys. Rev. Lett. 101, 

261104 (2008). 
[17] P. Blasi and P. Serpico, Pr eprint larXiT: 0904.0871 
[18] Fermi/LAT Coll., Preprint arXiv:0905.0025 



